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Abstract.  This  paper  describes  numerical  simulation  of  rarefied  nozzle  plume  impingements.  Two  different 
reservoir  pressures  400  kPa  and  4  kPa  axe  considered.  In  the  case  of  400  kPa,  the  simulation  of  the  nozzle 
flow  was  conducted  by  using  the  Navier-Stokes  equation,  and  then  the  analysis  of  the  plume  flow  was  carried 
out  by  the  DSMC  method,  employing  the  nozzle  exit  conditions  obtained  by  Navier-Stokes  equation.  On 
the  other  hand,  for  4  kPa,  both  the  nozzle  flow  and  the  plume  impingement  have  been  calculated  using  the 
DSMC  method.  Concerning  the  angle  between  the  nozzle  axis  and  the  flat  plate,  three  kinds  of  angle  are 
selected,  that  is,  90°,  45°  and  0°.  In  addition,  we  considered  the  case  where  there  exists  a  flat  plate  behind 
the  nozzle.  Simulated  results  have  been  compared  with  the  existing  experiments  for  the  pressure  and  shear 
stress  distributions  on  the  flat  plate.  A  good  agreement  between  the  DSMC  results  and  the  experiments 
are  shown.  In  the  case  of  the  oblique  and  parallel  impingements,  the  location  of  the  impingement  pressure 
peak  and  the  stagnation  point  shifted  upstream  with  increasing  rarefaction. 


INTRODUCTION 

Satellites  and  spacecrafts  are  usually  equipped  with  thrusters  for  attitude  control.  Exhaust  gas  from  the 
thruster  expands  into  the  space  vacuum  and  generates  a  huge  plume  in  the  space.  As  a  result,  the  plume 
may  impinge  on  the  satellite  body,  antennas  or  solar  panels.  This  impingement  causes  significant  problems 
of  surface  degradation,  contamination,  heat  loads,  disturbance  torque  and  so  forth.  Hence,  it  is  worth  while 
investigating  the  detailed  structure  of  the  plume  impingement  for  the  purpose  of  space  engineering. 

The  flat  plate  impingements  of  a  nozzle  plume  have  been  experimentally  examined  by  Legge  [1]  in  a  high 
vacuum  facility.  He  measured  pressure  and  shear  stress  on  a  flat  plate  for  normal  and  oblique  impingements, 
and  compared  with  the  impingement  model  [2]  based  on  Simon’s  plume  model  [3].  Concerning  the  numerical 
simulation,  the  direct  simulation  Monte  Carlo  (DSMC)  method  introduced  by  Bird  [4]  has  been  widely  used 
in  rarefied  gas  flows  with  the  development  of  computers.  The  problem  of  the  nozzle  plume  impingement  on 
a  perpendicular  surface  has  been  studied  by  using  the  DSMC  method  [5].  However,  there  are  no  numerical 
studies  concerning  parallel  and  oblique  impingements  on  a  flat  plate. 

Figure  1  illustrates  several  types  of  a  nozzle  plume  impingement,  where  (3  is  the  angle  between  the  nozzle 
axis  and  the  flat  plate  and  9um  the  maximum  turning  angle  of  a  stream  fine  at  the  nozzle  exit.  The  model 
field  means  the  area  where  Simon’s  plume  model  is  applicable.  The  model  field  is  dependent  on  the  angle  9um. 
Therefore,  there  exists  the  area  where  the  impingement  model  is  not  available,  especially  in  the  cases  of  (c)  and 
(d).  In  the  present  paper,  we  have  performed  the  numerical  simulation  of  nozzle  plume  impingements  on  a  flat 
plate,  and  calculated  the  impingement  pressure  in  the  area  where  the  impingement  model  is  not  applicable.  In 
the  case  of  the  normal  and  oblique  impingements,  we  also  compared  quantitatively  with  Legge’s  experimental 
data  for  the  pressure  and  the  shear  stress  distributions  on  the  flat  plate. 


NUMERICAL  METHOD 


A  plume  impingement  considered  in  the  present  study  is  such  that  gas  is  expanded  through  a  nozzle  into 
vacuum  and  then  interacts  with  a  flat  plate.  The  diameters  of  the  nozzle  throat  and  the  nozzle  exit  are  0.6  mm 
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and  4.75  mm,  respectively.  The  nozzle  length  is  7.75  mm  and  the  nozzle-divergent  angle  is  30°.  In  Figs.  2,  3 
and  4,  the  computational  fields  for  Fig.  1  (a),  (c)  and  (d)  are  shown,  respectively.  The  numerical  simulations 
for  the  reservoir  pressure  Po=400  kPa  (case  1)  and  4  kPa  (case  2)  are  conducted  in  the  present  paper.  For 
case  1,  the  continuum  theory  is  applicable  inside  the  nozzle.  Hence,  we  simulated  the  nozzle  flow  by  using  the 
Navier-Stokes  equation.  DSMC  simulation  of  the  plume  flow  is  carried  out  by  using  the  nozzle  exit  conditions 
obtained  by  Navier-Stokes  simulation.  On  the  other  hand,  for  case  2,  both  the  nozzle  flow  and  the  plume 
impingement  have  been  calculated  using  the  DSMC  method.  In  the  cases  of  Figs.  2  and  3,  the  distance 
between  the  flat  plate  and  the  center  of  the  nozzle  exit  L  is  set  to  40  mm  in  order  to  compare  the  simulation 
with  Legge’s  experiments.  Two  different  distances  of  L  are  selected  in  Fig.  4,  0  mm  and  15.5  mm. 


In  the  present  calculation,  nitrogen  was  adopted  as  a  test  gas  to  compare  with  Legge’s  experiments.  The 
flowfield  model  considered  here  is  such  that  molecules  impinging  on  a  flat  plate  suffer  diffuse  reflection,  and 
the  temperature  on  the  flat  plate  is  assumed  as  equal  to  laboratory  temperature  (300K).  Molecules  across  the 
symmetric  plane  suffer  specular  reflection.  Outflow  conditions  are  used  as  ambient  conditions.  Collisions  of 
molecules  are  simulated  using  the  variable  hard  sphere  (VHS)  model  [6].  The  Borgnakke-Larsen  statistical 
model  [7]  is  employed  for  the  calculation  of  the  energy  exchange  between  translational  and  rotational  modes, 
together  with  temperature-dependent  energy  exchange  probability  of  Boyd  [8].  The  computational  domain 
consists  of  9  to  14  subregions  for  efficient  computations,  with  total  cells  of  405,300.  The  total  number  of 
simulated  molecules  in  this  domain  is  9.0  x  105  in  the  case  of  the  oblique  impingement  (case  1). 
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FIGURE  1.  Impingement  types 


FIGURE  2.  Computational  field  for  Fig.  1  (a) 


FIGURE  3.  Computational  field  for  Fig.  1  (c) 


FIGURE  4.  Computational  field  for  Fig.  1  (d) 
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RESULTS  AND  DISCUSSION 


First,  the  numerical  simulation  of  a  nozzle  plume  is  conducted  without  considering  the  flat  plate.  Figures 
5  and  6  show  translational  temperature  and  number  density  contours  for  case  1  in  the  nozzle  and  the  vicinity 
of  the  nozzle  exit.  At  the  center  of  the  nozzle  exit,  Mach  number  is  5.8,  the  temperature  is  39  K  and  the 
Knudsen  number  is  1/2087.  It  is  difficult  to  make  the  cell  of  the  mean  free  path  length  because  the  vicinity  of 
nozzle  exit  is  still  a  continuum  flow.  In  the  present  DSMC  simulation,  therefore,  the  cell  which  is  larger  than 
the  mean  free  path  length  was  employed  near  the  nozzle  exit.  In  Figs.  7  and  8,  translational  temperature  and 
number  density  contours  for  case  2,  respectively,  are  illustrated.  The  nozzle  flow  simulation  for  case  2  indicates 
that  Mach  number  is  3.9,  translational  temperature  is  74  K  and  the  Knudsen  number  is  1/60  at  the  center  of 
the  nozzle  exit.  As  the  reservoir  pressure  is  decreased,  the  nozzle  flow  is  influenced  by  the  nozzle  wall  and  the 
boundary  layer  becomes  large.  Comparison  of  Fig.  6  with  Fig.  8  shows  that  the  backflow  of  the  nozzle  plume 
for  case  2  is  more  significant  than  that  for  case  1. 
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FIGURE  5.  Translational  temperature  contours  for  case 
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FIGURE  6.  Number  density  contours  for  case  1 
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FIGURE  7.  Translational  temperature  contours  for  case 
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FIGURE  8.  Number  density  contours  for  case  2 


Figure  9  shows  the  comparison  between  the  simulation  and  experiments  of  the  pressure  distribution  on  the 
flat  plate  for  the  normal  impingement  (/ 3  =  90°).  As  shown  in  this  figure,  there  is  an  excellent  agreement 
between  the  simulation  and  the  experiments  for  the  impingement  pressure  distributions.  It  is  indicated  that 
the  DSMC  method  is  valid  to  simulate  the  nozzle  plume  impingement  by  using  the  nozzle  exit  conditions  from 
the  Navier-Stokes  solution,  in  spite  of  large  cell  near  the  nozzle  exit.  In  the  present  paper,  the  numerical 
simulation  for  the  more  rarefied  conditions  than  case  1  is  conducted  to  investigate  the  effect  of  the  rarefaction, 
this  is  case  2.  Comparison  between  case  1  and  case  2  is  also  illustrated  in  Fig.  9.  Inside  the  nozzle,  for  case  2, 
the  flow  is  already  rarefied,  therefore,  Mach  number  at  the  nozzle  exit  for  case  2  is  smaller  than  that  for  case 
1.  As  a  result,  it  is  seen  that  the  impingement  pressure  for  case  2  is  smaller  than  that  for  case  1.  Figure  10 
shows  the  comparison  between  the  DSMC  results  and  Legge's  experiments  for  the  shear  stress  distributions 
on  the  flat  plate.  The  present  result  shows  a  good  agreement  with  Legge’s  experiments.  The  minimum  of  the 
impingement  shear  stress  (x/L  =  0)  means  the  location  of  the  stagnation  point.  For  case  2,  the  maximum  of 
the  calculated  results  is  slightly  smaller  than  that  of  case  1. 
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FIGURE  9.  Pressure  distributions  on  the  flat  plate  for  FIGURE  10.  Shear  stress  distributions  on  the  flat 
0  =  90°  plate  for  0  =  90° 


FIGURE  11.  Pressure  distributions  on  the  flat  plate  FIGURE  12.  Shear  stress  distributions  on  the  flat 
for  0  =  45°  plate  for  0  =  45° 


In  the  case  of  the  oblique  impingement,  the  angle  between  the  nozzle  axis  and  the  flat  plate  is  45°.  Figure  11 
illustrates  the  comparison  between  the  DSMC  results  and  the  experiments  for  the  impingement  pressure.  As 
shown  in  this  figure,  the  DSMC  result  for  case  1  indicates  the  feature  similar  to  the  experiments.  Comparison 
between  case  1  and  case  2  shows  the  location  of  pressure  peak  is  shifted  upstream  on  the  plate  with  increasing 
rarefaction  in  case  2.  Figure  12  shows  the  comparison  for  the  shear  stress.  The  DSMC  result  for  case  1  indicates 
a  good  agreement  with  the  experimental  data  at  x/L  >  0.45.  It  is  seen  that  the  location  of  the  stagnation 
point  is  between  x/L  =  0  and  the  location  of  the  corresponding  impingement  pressure  peak,  that  is  x/L  =  0.45 
for  case  1  and  x/L  =  0.05  for  case  2,  respectively.  Unlike  the  normal  impingement,  the  maximum  of  the 
impingement  shear  stress  become  large,  as  the  reservoir  pressure  is  decreased. 

Next,  we  consider  the  parallel  impingement.  Figure  13  shows  the  comparison  between  case  1  and  case  2  for 
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the  impingement  pressure  distributions.  Experimental  data  for  the  parallel  impingement  don’t  exist,  so  that 
the  comparison  of  the  simulation  with  the  experiment  is  not  shown  in  the  present  paper.  As  the  reservoir 
pressure  is  decreased,  the  maximum  of  the  impingement  pressure  shifts  upstream  on  the  plate.  That  means 
the  expansion  of  the  nozzle  plume  is  large  with  increasing  rarefaction.  In  the  case  of  the  oblique  and  normal 
impingement,  the  maximum  of  the  impingement  pressure  for  case  2  was  smaller  than  that  for  case  1.  In  the 
case  of  the  parallel  impingement,  on  the  other  hand,  the  maximum  of  the  pressure  for  case  2  is  larger  than 
that  for  case  1.  In  Fig.  14,  the  comparison  for  the  impingement  shear  stress  distribution  is  illustrated.  The 
shear  stresses  for  case  2  is  larger  than  those  for  case  1. 


FIGURE  13.  Pressure  distributions  on  the  flat  plate  FIGURE  14.  Shear  stress  distributions  on  the  flat 
for  0  =  0°  plate  for  0  =  0° 
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FIGURE  15.  Pressure  distributions  on  the  flat  plate 
for  0  =  —90°  (case  1) 
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FIGURE  16.  Pressure  distributions  on  the  flat  plate 
for  0  =  —90°  (case  2) 


Finally,  we  simulated  the  case  that  there  exists  the  flat  plate  behind  the  nozzle  (Fig.  4).  We  consider  L  =  0 
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mm  and  15.5  mm  as  the  distance  between  the  nozzle  axis  and  the  flat  plate.  Figure  15  shows  the  impingement 
pressure  on  the  flat  plate  for  case  1.  The  maximum  of  the  impingement  pressure  is  at  the  edge  of  the  nozzle 
(r  =  2.375  mm)  in  the  case  of  L  =  0  mm,  and  the  pressure  is  decreased  as  away  the  nozzle.  However,  the 
pressure  on  the  flat  plate  for  L  =  15.5  mm  is  almost  constant,  1.0  x  10-9.  There  doesn’t  exist  the  impingement 
pressure  at  0  mm  <  r  <  0.007  mm,  because  molecules  in  the  DSMC  simulation  don’t  exist  in  the  area  outside 
the  expansion  angle.  The  comparison  of  two  results  illustrated  that  the  impingement  pressure  for  L  =  0  mm 
and  15.5  mm  approach  the  value  P/Po  =  1.0  x  10-9  away  from  the  nozzle.  In  Fig.  16,  the  impingement  pressure 
for  case  2  is  illustrated.  The  comparison  of  Fig.  15  with  Fig.  16  shows  that  the  impingement  pressure  for  case 
2  is  larger  than  that  for  case  1,  because  the  number  of  molecules  scattering  behind  the  nozzle  is  increased, 
as  the  reservoir  pressure  is  decreased.  The  pressures  for  both  L  =  0  mm  and  15.5  mm  approach  the  value 
P/Pq  =  1.0  x  10“8,  away  from  the  nozzle. 


CONCLUSIONS 

The  DSMC  simulation  of  normal,  oblique  (/?  =  45°)  and  parallel  nozzle  plume  impingements  on  a  flat  plate 
yielded  the  following  conclusions:  By  combining  the  Navier-Stokes  equation  with  the  DSMC  method,  the 
numerical  simulation  of  the  nozzle  plume  impingement  was  conducted.  The  calculated  results  were  compared 
with  the  Legge’s  experiments  for  pressure  and  shear  stress  distributions  on  the  flat  plate,  and  a  good  agreement 
between  two  results  was  obtained.  It  may  be  mentioned,  therefore,  that  the  DSMC  method  is  valid  to  simulate 
the  nozzle  plume  impingement  by  using  the  nozzle  exit  conditions  from  the  Navier-Stokes  solution.  In  the  case 
of  the  parallel  impingement,  it  is  seen  that  the  normalized  impingement  pressure  and  shear  stress  become  large 
with  increasing  rarefaction.  In  addition,  the  simulation  of  the  case  that  there  exists  the  flat  plate  behind  the 
nozzle  (Fig.  4)  is  conducted.  It  was  revealed  that  the  impingement  pressure  approaches  a  value  away  from  the 
nozzle  independently  of  the  distance  between  the  nozzle  exit  and  the  flat  plate. 
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